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Abstract. We consider a spin-boson Hamiltonian which is generahzed such that the 
Hamihonians for the system (Hs) and the interaction with the environment (iJint) do 
not commute with each other. Considering a single-mode quantized field in exact 
resonance with the tunneling matrix element of the system, we obtain the time- 
evolution operator for our model. Using our time-evolution operator we calculate 
the time-dependent pointer states of the system and the environment (which are 
characterized by their ability not to entangle with each other) for the case that the 
environment initially is prepared in the coherent state. We show that our solution for 
the pointer states of the system and the environment is valid over a length of time 
which is proportional to n, the average number of bosons in the environment. We also 
obtain a closed form for the ofFdiagonal element of the reduced density matrix of the 
system and study the decoherence of the central system in our model. We show that for 
the case that the system initially is prepared in one of its pointer states, the offdiagonal 
element of the reduced density matrix of the system will be a sinusoidal function with 
a slow decaying envelope which is characterized by a decay time proportional to n; 
while it will experience a much faster decoherence, when the system initially is not 
prepared in one of its initial pointer states. 

Keywords: Spin-boson model, Jaynes-Cummings model, Decoherence, Pointer states of 
measurement. Foundation of quantum mechanics. 
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1. Introduction 

This is the second paper in the series of papers where we discuss pointer states of 
measurement and their evaluation. We refer to this paper as paper II in this series of 
papers. 

1.1. Foreword 

The pointer states of a system are defined as those states of the system characterized by 
their property of not entangling with states of another system [11 [21 13]. This condition 
is commonly referred to as the stability criterion for the selection of pointer states PQ |2] . 



Time-dependent pointer states of the generalized spin-boson model and etc. 



2 



In paper I |1] we discussed the pointer states of measurement and presented a general 
method for obtaining the pointer states of the system and the environment for a given 
total Hamiltonian defining the system-environment model. As we elaborately described 
in paper I, generally we should distinguish between the pointer states of a system and 
the preferred basis of measurement. We explicitly proved that the pointer states of a 
system generally are time-dependent and a preferred basis of measurement does not 
exist, unless under some specific conditions (discussed there in paper I) that the pointer 
states of measurement become time-independent. Moreover, time-dependent pointer 
states necessarily are not orthogonal amongst themselves at all times. Therefore, they 
cannot be considered as eigenstates of a Hermitian operator at all times. 

In paper I we also used our method in order to rederive the time-dependent pointer 
states of the system and the environment for the Jaynes-Cummings model (JCM) of 
quantum optics verifying the previous results for the JCM [SI El Cj. 

In this paper we study a spin-boson model [|] which is defined through the following 
total Hamiltonian 



This model basically is composed of a central spin-half particle (or other two-level 
system) surrounded by an environment of bosonic particles such as photons. In the 
Hamiltonian of the spin-boson model, represented by equation ([T]), we have considered 
an intrinsic tunneling contribution in the self-Hamiltonian of the system (proportional 
to the ax Pauli matrix), which can induce intrinsic transitions between the upper and 
lower states of the central system. Here Aq is the so-called tunneling matrix element. 
Also, it is assumed that the asymmetry energy in the self-Hamiltonian of the central 
system is negligible. Therefore, here we do not consider a contribution proportional 
to the (Jz Pauli matrix in the self-Hamiltonian of the system [§j The second term in 
equation ([T]) represents the self-Hamiltonian of the electromagnetic field; where we have 
considered a single-mode quantized field with the frequency of u for the environment. 
Also the third term, with g as the spin- field coupling constant, represents the interaction 
between the central spin-half particle and a single-mode quantized field; which in fact 
is the quantized form of the famous — ^t.B Hamiltonian due to the interaction between 
a particle of magnetic dipole-moment fi and a magnetic field B. 

In the Hamiltonian of equation ([T]), if we switch the ax and az operators and 
apply the rotating-wave approximation, we would obtain a Hamiltonian which would 
precisely look like the Hamiltonian of the Jaynes-Cummings model (JCM). However, the 
Hamiltonian that would be obtained in this way, will not describe the Jaynes-Cummings 
model of quantum optics. The above point is because of the fact that when considering 

I For a serious review and analysis of spin-boson models in different regimes the interested reader can 
refer to the seminal article by Leggett et al. |8] or the book by Weiss [9] . Also a brief while very useful 
review of the model can be found in chapter 5 of Schlosshauer's book [T]. 

§ However, one can easily verify that modifications due to such contribution can be done quite easily, 
as such a term commutes with the interaction between the system and the environment, represented 
by the third term in equation (fll). 




(1) 
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our model, represented by the Hamiltonian of equation ([T]), a very special meaning is 
attributed to the eigenstates of the operator; they are the upper and lower states 
of the two-level system. In fact, this attribution of the upper and lower states of the 
two-level system to the eigenstates of the operator will be considered everywhere in 
our calculations; as is also considered in calculations for the JCM [101 [H]. Therefore, as 
we will see in the following section, when considering the rotating-wave-approximation 
with the exact resonance condition the Hamiltonian of equation ([T]) will lead to a time- 
evolution operator that is very different from the one that we know for the JCM. As 
a result, as we will see, there will be clear differences between the physics which arises 
from our spin-boson model and the one that we know from the JCM. 

In the paper by Leggett et al. [8] they considered a general form of the Hamiltonian 
given by equation ([T|, where the environment can be represented by a spectral density 
J{u}) (rather than considering a single-mode quantized field). They used the "influence- 
functional" method of Feynman and Vernon |T2] to obtain general expressions for 
P(t) = (azit)) in the form of a power series in A. However, the general expressions 
they obtained for P{t) in terms of the spectral density function J(w) were exceedingly 
cumbersome to calculate in most regimes. So, they assumed a nonresonance regime, 
which was quite different from the resonance regime that we will consider here in this 
article 0. 

The model represented by the Hamiltonian of equation ([T]) can also be studied in 
the framework of the Born-Markov approximation in order to obtain an approximate 
master equation for the evolution of the reduced density matrix of the system [9]. The 
master equations obtained in this way are valid only in certain regimes; and moreover, 
one often may need to resort to numerical computation in order to be able to solve them. 
However, the main purpose of this paper is (1) to obtain the time- dependent pointer 
states of the system and the environment, as well as expressions for the evolution of 
the reduced density matrix of the system in the exact resonance regime and for an 
environment initially prepared in the coherent state; and (2) to obtain approximate 
expressions in closed form for the evolution of the off-diagonal elements of the reduced 
density matrix (for the case that the environment initially is prepared in a coherent 
state with a large average number of photons), which can be used in order to study the 
decoherence of the central system in an analytical way. 

1.2. Review of our method for calculation of pointer states 

Our goal is to find the pointer states of the system, as well as their corresponding 
pointer states from the environment for an arbitrary total Hamiltonian defining the 
system-environment model. In order to find pointer states, we look for system states 
\si{t)) such that the composite system-environment state, when starting from a product 
state |sj(to))|-E'o) at t = 0, remains in the product form \si{t))\Ei{t)) at all subsequent 
times t > 0; i.e. pointer states must satisfy the aforementioned stability criterion. 
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Now consider a two-state system S with two arbitrary basis states |a) and \b), 
initially prepared in the state 

\^^{to)) = a\a) + I3\b) with + |/3|2 = i; (2) 

and an environment initially prepared in the state 

oo 

\<^'{to)) = J2^nM, (3) 
n=0 

where {|(/3ri)}'s are a complete set of basis states for the environment. For the two-state 
system with the two basis states \a) and |6) we can consider a time evolution operator 
for the global state of the system and the environment given by 

f/tot(t) = S,\a){a\+S2\a){b\+Ss\b){a\+£,\b){b\; (4) 

where the £^j's are some generally time- dependent operators in the Hilbert space of 
the environment, and depend on the total Hamiltonian defining the system-environment 
model. 

Using equations (2) to (4) we can write the global state of the system and the 
environment as 

oo 

\^'°\t)) = [/tot(t). ia\a) + m) ® (E ^nM) 

n=0 

= A(t) |a) + B(t) (5) 

oo 

with A{t) = E Cn{a£i{t) + M 

n=0 

oo 

and Bit) = c„{a4(t) + /3^4(t)} 

n=0 

Now, for the global state of the system and the environment, given by equation (|5]), we 
observe that if for some initial states of the system and the environment the two vectors 
A(t) and B(t) of the Hilbert space of the environment turn out to be parallel with each 
other, i.e. if 

A{t) = G(t)B(t), (6) 

with G{t) as a scalar which generally may depend on time, then those initial states 
of the system and the environment will not entangle with each other, and hence they 
can represent the initial pointer states of the system and the environment. In fact, if 
for some initial states of the system and the environment the condition represented by 
equation ^ is satisfied, the global state of the system and the environment (equation 
(Isj)) can be written in a product from as 



= {G{t)\a) + \b)} X Cn{aSs + PS,} (7) 

n=0 

in which case pointer states can be realized for the system and the environment given 
by 

I ± (t)) = Af {G{t)\a) + \b)} and (8) 

oo 
n=0 
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In the above equation we have represented the pointer states of the system by | ± (t)) 
and those of the environment by |$-|-(t)). Also, JV is the normahzation factor for the 
pointer states of the system. 

In general, an arbitrary choice of c„'s will not necessarily yield an A(t),B(t) pair 
that remain parallel for any choice of a and f3; therefore an arbitrary choice of c„'s does 
not necessarily correspond to a pointer state for the environment. Nonetheless, for a 
given set of c„'s there may exist two sets of complex numbers for a and /3 with two 
corresponding values for the scalar function G(t) such that equation ^ is satisfied. In 
the following example, we show that for an assumed set of c„'s, values of a and (3 exist 
such that equation (|6| is approximately satisfied in the limit of large number of degrees 
of freedom for the environment fill 



2. Calculation of the time-evolution operator 



In order to calculate the time-evolution operator in the interaction picture for the 
Hamiltonian given by equation ([T]), first we need to have the corresponding Hamiltonian 
in the interaction picture. It can easily be calculated as 

i^int(t) = {a. cos(Aot) - ay sin(Aot)}{gate*-* + g*de-'^'}. (9) 

Here, the commutator of Hmtit) and Hintit' 7^ t), i.e. [Hintit) , Hmtit' 7^ t)] with Hmtit) 
given by equation (|9]), is not a function of a constant number. This in fact can make 
the evaluation of the time-evolution operator quite difficult [13] . 

In parallel with paper I we consider the general form given by equation Q for the 
evolution operator of the global spin-field system. For such time-evolution operator in 
the interaction picture, which satisfies the Schrodinger equation ih-^Uit) = HmtUit), 
we have 

I j = Hmtit) (I I) = {ga^e^-* + g^ae"^-*} (10) 

/ £1 cos(Aot) -|- iSs sin(Aot) S2 cos(Aot) + 184 sin(Aot) 
\ —£3 cos(Aot) — i£i sin(Aot) —£4 cos(Aot) — i£2 sin(Aot) 
Now, we assume the transition matrix element Aq to be in resonance with the cavity 
eigenmode u and we use the rotating-wave approximation (RWA) [TUl E] (just as is 

II In fact, by looking at equation ^ we notice that in practice we can expect some states of the system 
and the environment to keep their individuahty and not to entangle with each other even if they can 
satisfy our condition (given by equation Q) only in a fraction of the Hilbert space of the environment 
where the c„ coefficients are not negligible. This of course will involve assuming some approximations 
in obtaining pointer states. However, as we will show in this paper, in the end we can define a measure 
for the degree of entanglement between the states of the system and the environment, which after its 
calculation for the pointer states which we obtain for our model we can know exactly in which regimes 
our pointer states are valid and will not entangle with the states of another system. For example, this 
way we will show in this paper that the pointer states which we will obtain for our spin-boson model for 
an environment initially prepared in the coherent state are valid (i.e. will not entangle with the states 
of any other subsystem throughout their evolution with time) up to times of the order n/g; where n is 
the average number of photons in the coherent state of the environment. 
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assumed in the conventional Jaynes-Cummings model of quantum optics [ID])- So, by 
entering the resonance condition Aq = u and using the rotating wave approximation (i.e. 
disregarding the higher-frequency terms which contain e±«(^+^o)i^ ^j-^g above equation 
will simplify to the following set of four equations 



iSi = 


fa. 


-4) 


+ 




182 = 




-4) 


+ 




iSs = 


fa. 


-4) 






i£i = 


fiS2 


-£4) 







(11) 



In order to solve the above set of coupled differential equations, we proceed as 
follows. First, we take a derivative with respect to time of the first equation. By 
replacing Si and S3 from the first and the third equations in the resulting equation we 
find 

tii = ^^{il + 2N)Si-S3}, (12) 

where N = a^d is the number operator. Similarly, by doing the same procedure on the 
third equation for S3 we find 

-{il + 2N)S3-Si}. 



iSs 

Next, we define S. 



2 

and 

S+-\- = Si + Ss 



as follows 
and 



S-\— — Si — S3. 



(13) 



(14) 



By adding and subtracting equations (12) and (13) we find 

-|gpiV^++ and 
-\g\\N + l)S^_ . 



S 



++ 



(15) 



These equations for S+^ and can simply be solved to find 

S++ = sin(|g|\/i^ t)A + cos(|g|\/iV t)B and 

= sin(|g|\/iV + 1 t)A' + cos(|g|\/iV + 1 t)B', (16) 

where A, A', B and B' are some time-independent operators (rather than constant 
numbers)., which will be found from our initial conditions in the following paragraphs. 
Here we note that since these coefficients generally are some time-independent operators 
rather than constant numbers, and they do not necessarily commute with the number 
operator iV, we must have them on the right-hand side of the sin and cos functions 
(rather than having them on the left-hand side). Only in this way the solutions in 
equation (16) will satisfy equation (15). 
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Now using equations ( 14 ) and ( 16 ) , we can obtain the operators £i and £3 as follows: 



Si = -{sindglVA^ t)A + cos(|g|V t)B 



+ sin(|g|ViV + 1 t)A' + cos(|g|ViV+ 1 t)B'} and 
1 



£, = ^{sindgl ViV t)A + cos(|g|ViV t)B 



(17) 



-sin(|g|VA^ + 1 t)A' - cos(|g|ViV+ 1 t)B'}. 
In quite the same manner we can calculate S2 and £^4 as follows 

£2 = ^{sin(|g|\/iV t)C + cos(|g|\/iV t)D 



+ sin(|g|VA^+ 1 t)C" + cos(|g|VA^+ 1 t)D'} 
4 = ^{sin(|g|\/i^ t)C + cos(|g|\/iV t)D 



and 



-sin(|g|VA^ + 1 t)C'-cos(|g|ViV+ 1 t)D'}; 



(19) 



(20) 



where C, C, D and D' also, generally are some time-independent operators, which will 
be determined from our original set of equations (11 ) and the initial conditions on {Si}^s. 

In order to obtain the eight operator coefficients which appear in our expressions 
for {^j}'s, first we note that the time-evolution operator, given by equation (4|, must 
satisfy the initial condition U(t = 0) = where Is = \o){a\ + \b){h\ represents the 

identity operator in the Hilbert space of the system and Is is the identity operator in 
the Hilbert space of the environment. This means that we must have 

^i(O) = ^4(0) = is and ^2(0) = 4(0) = (21) 

From the above initial conditions and equations ( [l7) ) to (21) we easily find four of the 
coefficients as follows 



B = B' = D = le 



and 



D' = -If 



(22) 



In order to find A and A' we proceed as follows. First, we use equation (11 -a) to 
obtain £3 as follows 

(g*a - go^% = 2iti - (g*a + ga^% . (23) 



Replacing £1 and £1 from equation (17) into the above equation, it reads 
(g*a - ga))£z = i\g\\/^ {cos(|g|\/i^ t)A - sin(|g|\/i^ t)} 



+i|g| VA^ + 1 {cos(|g| ViV + 1 t)A' - sin(|g| ViV + 1 t)} 
g*a -I- go) 



(24) 



){sin(|g|V t)A + cos(|g| V t) 



+ sin(|g|VA^ + 1 t)A' + cos(|g|VA^+ 1 t)}. 
At t = the above equation reduces to 



2|g|VA^ A + i\g\\JN + 1 A' -{g*a + ga 



it) 4(t = 0) = 0.(25) 
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A \n) +i|g|VA^ + 1 A' \n) = g'y/n \n- 1) + g\/n + 1 |n+ 1).(26) 
Next, we use equation (11 -c) to obtain £i as follows 

(gat - g*a% = 2z4 + {g*a + ga^% . (27) 
Replacing E-^ and from equation (19) into the above equation, it reads 

(gat - g*a)^i = i\g\\[f^ {cos(|g|\/i^ t)A - sin(|g|\/i^ t)} 

+i\g\^N + 1 {-cos(|g|\/iV + l t)i' + sin(|g|\/iV + l t)} (28) 
o-*a + gat 



-){sin(|g| V t)A + cos(|g| V t) 



+ ( 

-sin(|g|\/A^+ 1 t)A' - cos(|g|\/A^ + 1 t)}. 
At t = the above equation reduces to 

i\g\\[^ A - i\g\^N + 1 A' = (gat - g*a) Si{t = 0) = (gat - g*a) 
Operating this last equation on \n) we have 



(29) 



N A\n)- i\g\\/N + 1 A' \n) = g^n + 1 |n + 1) - g*y/n \n - 1).(30) 



Finally, we use equations (26) and (30) to obtain the coefficients A and A'. Assuming 

(31) 



g to be real and then adding equations (|26|) and (|30|) we find 

A\n) = - 



ln + 1 



N 



Comparing the above equation to 



n + 1) = —i\n -\- 1) . 

at \n) = —i\n-\- 1) we find A as 



A 



—I 



at. 



Similarly, by subtracting equation (30) from equation (26) to find 

A' \n) = —i. / — \n — 1) = —i\n — 1) 

\ N+1 

and comparing the above equation to — a In) = —i In — 1) we find A' as 



A' 



a . 



(32) 



(33) 



(34) 



'A^ + 1 

One should pay attention to the order that the operators appear in equations (32) and 
(34); as they are not commuting operators. 



By doing exactly the same procedure on equations (11-b) and (11 -d) we would find 
the operator coefficients C and C' as follows 



C 



-A 



and 



C 



a = A'. 



N 



(35) 
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Now that we found all the operator coefficients, we can replace them back in equations 
(17) to (20) and write the {Si} in their final form 

1 



Si = ^{-2sin(|g|VA^ + cos(|g|\/A^ t) 

1 



isin(|g|vA^ + 1 t) 



a + cos(|g|VA^ + 1 t)} 



'N + 1 

4 = ^{-«sin(|g|\/i^ t)^ a"^ + cos(|g|VA^ t) 

yN 
1 



+isin(|g|VA^ + 1 t) 



N + 1 



a-cos(|g|VA^ + 1 t)} 



S2 = ^{isin(|g|\/i^ + cos(|g|\/i^ t) 

\ N 



isin(|g|ViV + 1 t) 



N + 1 



a - cos(|g| yiV + 1 t)} and 



4 = -{isindglViV t)^= + cos(|g|\/i^ t) 

1 



+zsin(|g|ViV + 1 t) 



a + cos(|g|VA^ + 1 t)} 



(36) 



(37) 



(38) 



(39) 



'N + 1 

One can verify that the above set of operators satisfy the unitarity of the time-evolution 
operator = UU^ = 1. 



3. Calculation of the time-dependent pointer states of the system and the 
environment 

In order to obtain the time-dependent pointer states of the system and those of its 
environment for our model, represented by the Hamiltonian of equation ([T]), we assume 
the field to be initially prepared in the coherent state |z/) 

|$fieid(to)) = W) = Z^'^^l^)'^ ■wi^'^ = — ' (40) 

n=o Vnl 

where |z/p = fi is the average number of photons in the coherent state, and u = \v\e~'^'^. 
In this section we will show that in the regime that we are considering (i.e. the exact- 
resonance with the RWA regime) and for the environment initially prepared in the 
coherent state, in the limit of a large average number of photons n — ?■ 00 we must have 
pointer states for the system (the central spin) given by 

I + (t)) = -zcos(f + -^)|a) + sin(f + -^)|6) and 
2 4v?T, / 4v^ 



- it)) = zsin(f - -^)|a) + cos(f - -^)|6), (41) 
2 4v^ 2 4v^ 
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where \a) and \b) are eigenstates of the Pauh matrix. 

We make the usual assumption that there exists no correlations between the system 
and the environment at t = 0. So, we consider the following initial state for the total 
composite system 

\ij'°\to)) = {a\a) + f3\b))(^T.Cn\n) with \a\^ + \f3\^ = 1. (42) 

n=0 

For pointer states the two vectors A{t) and B(t) of the Hilbert space of the environment 
must be parallel with each other. Therefore, for pointer states due to our condition, 
given by equation ([6]), we must have: 

^c4aA + /34} M = Git) xJ2cn{a£3 + ^£4} M- (43) 

n n 

For the {Si}, given by equations (36) through (39), we can write 

£M) = Mn)\n + 1) + Mn)\n) + fM\n - 1), (44) 
where fij's (with z = 1, 2, 3, 4 and j = 1, 2, 3) are given by 

fii{n) = f3i{n) = -f2i{n) = -Ui{n) = y sin(gtyri^TI) = fi{n), 

fi2{n) = /42(n) = ^(cos(gtv^) + cos{gtVn + 1)) = /2(n), 

/22(n) = /32(n) = ^(cos(gtv^) - cos(gtVn+ 1)) = fsin) and (45) 

— t 

fu{n) = f23{n) = -fssin) = -hsin) = y sin(gtv/^) = f[{n). 



Using equation (44) our condition, given by equation (43), becomes 

00 

E Cn{aifi^{n)\n + 1) + /i2(n)|n) + Mn)\n - 1)) 

n=0 

+f3{f2i{n)\n + 1) + f22{n)\n) + f2s{n)\n - 1))} 

00 

= Git) X E Cn{a{f3i{n)\n + 1) + f32in)\n) + fMln - 1)) (46) 

n=0 

+/3(/4i(n)|n + 1) + /42(n)|n) + Mn)\n - 1))}. 

Now, since the number states {|^)} are a complete set of basis states for the 
environment, for the initial pointer states we can open the summations in equation 



(46) and equalize terms from the two sides of this equation which correspond to the 



same number state \n). After using equation (45), equation (|46|) would simplify as 



^ _ (a - /3)cnfi{n) + acn+if2{n + 1) + (icn+ih{n + !) + (« + P)cn+2f'i{n + 2) _ 
(a - /3)c„/i(n) + acn+ifsin + 1) + /5c„+i/2(n + !)-(« + /3)c„+2/((n + 2) ' 

for all n. (47) 

The above result for G(t), which generally depends on n, would contradict our initial 
assumption of the two vectors A(t) and B(t) being parallel to each other, unless we can 
find certain initial states for the system for which G{t) turns out to become independent 
of n; since as we discussed, for pointer states all components of the vector A {A^ s) must 
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be related to their corresponding components from B s) through the same scalar 
factor G (see equation (6)). So now we should seek for those particular initial states 
of the system which can make G{t) independent of the index n of the states of the 
environment. 



For an initial coherent field (equation (40)) we have: 



Cn+i = CnC ^'^\/-TT and Cn+2 = c„e 2»y — n Moreovcr, in the limit of a large 

^ V "+1 "-^^ " ^(n+l)(n+2) ' ^ 



average number of photons n — )■ oo we can replace the factors y and ^J by unit}- ^ 
since for n — )■ cxd the Poisson distribution of the coherent field is extremely sharp (wim 
n at the center) and hence, for n — )■ oo and n ^ n we have ^/^Jj ~ 1 and ^ 1, 

while for n being far from n the Cn coefficient is negligible. So, the corresponding terms 
(of n being far from n) do not have any contribution in the summations of equation 



(46). As a result, equation (47) for G{t) can be further simplified to 



^ (a - /3)/i(n) + ae-^^Mn + 1) + /3e--^/3(n + I) + (a + P)e-^'^ f[{n + 2) 
(a-/S)/i(n) + ae-^^/3(n + l) + /3e-^^/2(n + l)-(« + /3)e-2i^/{(n + 2)' ^ 

Replacing the functions from equation (lisj), the above equation reads 



G{t) = {(« + (5)e-''P cos{gtVn + 1) - z(a - /?) sin(gtVn + 1) 
+ {a - /3)e-^'^ cos(gtV^ + 2) - i{a + (3)e-^''^ sin(gtV« + 2)} 
^ {(a + /3)e-^'^ cos(gtVn + 1) - i{a - (3) sin(gtVn + 1) 
+ (/3 - a)e~''^ cos(gtV^ + 2) + i(a + (3)e-'^''^ sin(gtVri + 2)}. (49) 

In order to obtain the pointer states of the system, we should look for the probable 



initial states of the system (represented by the coefficients a and /3 in equation (42)) 
which can make the expression (in the above equation) for G{t) independent of the 
index n of the states of the environment. On the other hand, by looking at equation 



(49) we realize that ii a — (3 = ±(a + /3)e~*'^ the expression for G{t) will be considerably 
simplified. In what follows we show that for a — /3 = ±(q; + /3)e~*'^, which is equivalent 
to the initial conditions for the state of the system given by 

= — i cos((y9/2) and /3+ = sin((y9/2) (for the plus sign) or 
q;_ = 'isin((y9/2) and /3_ = cos(v9/2) (for the minus sign), (50) 



G{t) of equation (49) will be independent of the states of the environment; provided 



we have a large average number of photons in the field n — )■ oo. Therefore, the initial 



conditions of equation (50) correspond to the initial states of the system which do not 
entangle with the states of the environment. After that, we will obtain the time evolution 
of these initial pointer states and followed by that we obtain the corresponding pointer 
states of the environment. 

^ This approximation has been used by Gea-Banacloche and other people [51 [5] in the study of the 
Jaynes-Cummings model of quantum optics. In fact, for the Jaynes-Cummings model it has been shown 
that an average number of photons only as large as twenty is enough to make this assumption a good 
approximation [5]. 
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For a — j3 = ±{a + j3)e the expression in equation (49) for G{t) simphfies to 



G(t) = (51) 

The above expression can be written as 



Git) 



|g-iV./2 gT^(v^^+V^) j|gi</p/2 g±^(V^-v^^) ipg-i</p/2 g± ^ (v^^i+T-v^^) j 
|gjip/2 g±^(V^^-V^I+T) -p g-j(^/2 g±^(V^^+T-V^I+2)|' 



(52) 



From the Taylor series expansion of -y/n + 2 — + 1 about n (the average number of 
photons in the environment), given by 

/ / 1 (n -\- 1 — n) , ^ 

^/^^-^/^^=-^-^———^ + ..., 53 



we notice that in the hmit of a very large average number of photons we can replace 
A/rzr+^ — i/n + 1 by [5] ; since for very large fi all terms containing the index n, 
which appear after the first term, are negligible; and the series is convergent. So, in the 



classical limit of n — )■ oo we can rewrite equation (52) for G{t) as 

|gi(p/2 g±j5l -|- g-^'/'/S g^J^l} 



G(t) 



|gjl/j/2 g='=4y5i ip g-i'/'/S g^4v1i} 



— icot(— H — ^7=) for the first sign 
2 AJn 



or =itan(— ^^) for the second sign; (54) 

2 4vn 



which clearly is independent of the index n of the states of the environment. 

In appendix A by calculating the degree of entanglement between the states of the 
system and the environment for the pointer states which will be obtained from the above 
result, we will show that this result is valid over a length of time which is proportional 
to n, the average number of photons in the field. 



The result of equation (54)) simply means that for the initial states of the system 
given by 

I + {to)) = -i cos((^/2) |a) + sin((^/2) \h) and | - (to)) = i sin((^/2) \a) + cos((^/2) |6) (55) 

the states of the system and the environment will not entangle with each other. 
Moreover, using equation ([8| which gives us the general time evolution of the pointer 



states of the system; and G{t) of equation (54) (which is independent of the index n of 



the states of the environment) we can find the time evolution of the pointer states of 
the system as follows 

| + (t))=Ar+{-.cot(^ + -^)|a) + |6)} and 

|-(t))=Ar_{ztan(f --^)|a) + |6)}; (56) 
2 A^n 
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where Af^ and A/"- are the normahzation factors for the \ + {t)) and | — (t)) states 
respectively. It is easy to verify that 



7V+ = sin(^+(t)) and A/l = cos(^_(t)), where e±{t) 



So, we can rewrite equation ( |56| ) as 

I + (t)) 

I - (t)) 



^cos(- + -^)|a) +sin(- + —^)\b) 



1 ± _gL 

2 Ax/fi' 



and 



(57) 



■ ■ t"^ g^ M \ , t"^ g^ MM 

2 4v^ 2 4y/n 



(5J 



which is the same as equation (41); Q.E.D. 

Next, we obtain the corresponding pointer states of the environment. Using 
equations ([s]) and (44) we have 

oo 
n=0 



= Af^' Cn{«± [f3i{n)\n + 1) + f32{n)\n) + fUn)\n - 1)] (59) 

n=0 

+/3± [/4i(n)|n + 1) + /42(n)|n) + /43(n)|n - 1)]}, 

where in the above equation a± and P± are those of the initial pointer states of the 
system given by equation (50). Let us first obtain |$_|_(t)); i.e. the pointer state of 
the environment corresponding to the | + (t)) state. Replacing the fij functions from 
equation (45), a+ and /3+ from equation (50) and A/V from equation (57) we have 



-{- sin(gtVn + 1) e"*'^/^ \n + 1) 



-2sin(^+(t)) 
+ Hcos(gtv^) e"^/^ + icos{gt^yn + l) e~"^/^ ] \n) 

+ sin(gtv^) e*'^/^ \n - 1)}. (60) 

Using c„±i ~ Cne^**^ for the coherent field and in the limit of n — )■ oo, the above relation 
can be written as 

oo 

n=o 2zsm(fc^+(tjj 
+ [cos(gtv^) e*^/2 - cos(gty^5n) e"^^/^ ] 
+zsin(gtVri + 1) e"''^/^ } \n) 

„^o2.sin(Mt)) ^ 

which can easily be simplified into the following final result for |$+(t)) 

oo 



|?7,). 



(61) 



(62) 



ra=0 



Following exactly the same procedure one can also find |$_(t)) as follows 

oo 

|$_(t)) = Y,Cn e'f^v^+v^) \n). 



(63) 



n=0 
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We should also mention that the pointer states of the system at t = to (equation 
(55)) are orthonormal and hence, they form a complete basis set for the state of 
the system. Therefore, the evolution of any initial state of the system \ipsito)) = 
a' I + (to)) + I — (to)) in contact with an initial coherent field |z/) can be expressed as 
a linear combination of the evolution of | + {to))\i>) and | — {tQ))\i'), and in the following 
diagonal form: 

{a' I + (to)) + I - (to))) W) ^ a' I + it)) |<|.+(t)) + I - (t)) |<l>_(t)). (64) 



4. State preparation at specific times 



One interesting feature of the pointer states of the system, given by equation (58), is 
that at specific times they coincide with each other. In fact, by looking at equation (56) 
we notice that at those times for which tan(| — j^) = — cot(| + ^^) the | ± (t)) states 
are equal to each other. One can easily verify that at ti = [An + l)7rv^/g (with n = 
0, 1, 2, ...) both of the | ± (t)) states will be in the common state given by 

I ± (tO) = ^sin(| - ^)|a) + cos(| - ^)|6); (65) 

while at t2 = (4n — l)'n-y/^/g (with n = 1,2, ...) the | ± (t)) states will be in the common 
state given by 

I ± (t,)) = 2sin(| + ^)|a) + cos(| + ^)\b). (66) 

The state preparation at these specific times basically means that whatever is the 
initial state of the system, at these specific times the states of the system and the 
environment are not entangled to each other and the system can be represented by a 



well-defined state of its own (see equation (64)). Moreover, as we see, these specific 
states clearly depend on the phase ip of the initial state of the coherent field. The 
same kind of phenomenon was also discovered in the simpler Jaynes-Cummings model 
of quantum optics by Gea-Banacloche [3] in 1991. 



5. Consequences regarding the decoherence of the central spin 

In this section we will use the pointer states of the system and the environment, which 
we already obtained in section 3, in order to obtain a closed form for the coherences of 
the reduced density matrix of the system. However, prior to that let us use the time- 
evolution operator, which we already obtained in section 2, to calculate the offdiagonal 
element of the reduced density matrix of the system in a more precise form. 



5.1. General expressions for the evolution of the state of the total composite system 
and the reduced density matrix of the system 
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Using equations ([s]) and (44) to obtain \i'tot{t)), we can write 

oo 

l^tot(^)) = c„{(a/ii(n) + f3f2i{n))\a,n + 1) + (a/i2(ri) + I3f22{n))\a,n) 

3 

+ (a/i3H + /3/23H)|«, ^ - 1) + («/3iH + Phi{n))\b, n + 1) (67) 
+ («/32H +/3/42H)|fe,n) + (a/33(n) +/3/43(n))|6,n- 1)}, 



n=0 



where /j/s are given by equation (45). Replacing these functions from equation (45) 
the above relation can be simplified into the following general form 

oo 

\Aot(t)) = 'Y{ca,n(t)\a,n) + Cb,n(t)\b,n)), with 



n=0 



Ca,n(t) 

a — P 



—i a -\- (3 

— (a - /3)c„_i sin(gtA/n) + c„[( — ^ — ) cos(gtyn) 



+ ( — — ) cos(gtV^ + 1)] - - (a + /3)c„+i sin(gtVn + 1) 



(68) 



Cfe,n(t) = 

a — /3 



— i a -\- (3 

y (a - /3)c„_i sin(gtyn) + c„[( — ^ — ) cos(gtA/n) 



, 2 ) cos(gtA/?^ + 1)] + - (a + /3)c„+i sin(gtV^ + 1). 

We can do the trace operation over the basis states of the environment, and obtain 
the reduced density matrix of the system S as follows 

oo oo 

p'it) = j:{n\p'°\t)\n) = E(^lV^tot(t))(^tot(t)|n) 



n=0 



n=0 



E ( |Ca,n(t)r \a){a\ + Mt)\^ \b){b\+Ca,n{t)cln{t) \a){b\+C.C. ). (69) 



n=0 



So, in the basis of the eigenstates of (i.e. in the basis of the \a) and \b) states) the 
elements of the reduced density matrix of the system must be given by 



Puit) = E Ca,n{t) Cb,ni^) = Ca,0 o + Ca,l i + Ca,2 4 2 + ... and 
n=0 

oo 

PnW = l-pf2W = E|Ca,nWr (70) 



n=0 



Replacing Ca,n{t) and Cb^n{t) from equation (68) in the above equation, after some algebra 
one finds 



pUt) = 7/o(t) + + \f2{t) + \*Mt) + (A - x*)m 

Puit) = igoit) + Sg^it) + \g2it) + \*gl{t)- 
where in the above equations the coefficients 7, 5 and A are given by 



and 



(71) 



7 



and 



5 



|a + /3p 



and 



A = -(lap - + aP* - Pa*). 



(72) 
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Also the fi 

oo 

fo{t) = E 



n=0 



m = E 

n=0 



f2{t) = E 

n=0 



Mt) = E 

n=0 



/4(t) = E 

n=0 

oo 

9o{t) = E 

n=0 



9i{t) = E 

n=0 



^72(t) = E 

n=0 



t) and gi{t) functions are given by 

|c„_i|^sin^(gtv^) + z[c„_ic* + c*_iC„] 

X sm{gty/n) cos{gt\/n + 1) — |c„p cos^(gtv^?2 + 1) ), 
|c„|^ cos^(gtv^) - i[c„c*_^i + c*c„+i] 

X cos(gtA/n) sin(gtV^ + 1) - |c„+ip sin^(gtA/^ + 1) ), 
-ic„_ic; sin(gty^) cos(gtyn) - c„_ic;_^_i 

X sin(gtA/n) sin(gtA/^ + 1) - iCnC*^+i sm{gt\/n + 1) cos(gtA/nTT)), 
zc„c*_i sin(gtyn) cos{gt^/n) + c„+ic*_i 

X sin(gi(:v^) sm(gt\/n + 1) + icn+ic^ sm(gt\/n + 1) cos (gtVn + 1)), 
|c„p cos(gi-\/n) cos(gtA/n + 1)), 

|c„_i|^sin^(gtyri) +2[c„c*_i 

-c*c„_i] sin(gtA/n) cos(gtA/^ + 1) + |c„|^ cos^(gtV?T- + 1) ), 
|cn|^ cos^(gtv^) + i[cncl_^i - c*c„+i] X cos{gty/n) sm{gtVn + 1) 

+ |c„+i|^ sin^(g)f:A/n + 1)) and 
-zc„_ic; sm{gty/n) cos{gt^/n) + c„_iC*+i 



X sm{gty/n) sin(gtv^rH-T) + |c„|^ cos(gtA/n) cos(gtv^ n + 1] 
+2C„c*_|_]^ sin(g)f:A/ri + 1) cos {gtVn + 1)) 



(73) 



In Figure 1 we have used equation ( 71 ) to plot the evolution of the population 
inversion W{t) = Piiif) — p22{t) for the case that the two-level system initially is prepared 
in the upper level. Here the evolution of the population inversion is characterized by 
an oscillating envelope, which dominates the fine oscillations which show themselves as 
collapses and revivals of the atomic inversion in the simpler example of the JCM. 

Now, let us obtain the coherences of the reduced density matrix of the system in 
another way by using the pointer states of the system and the environment which we 
obtained in section 3. As we will see, in this way not only we can obtain a closed form 
for the coherences of the reduced density matrix of the system, but also we can acquire a 
better understanding regarding the characteristics of decoherence of the central system 
in our model. 



For \il)tot{t)) given by equation (64) the reduced density matrix of the system f)s{t) 
can be calculated by tracing over the environmental degrees of freedom to obtain 



Psit) = laf X I + (t))(+(t)| + X I - (t))(-(t)| + a'P^ 
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x| + X + /3V* X I - {t)){+{t)\ X ($+(t)|$_(t)). (74) 

So, in an arbitrary basis \a) and \b) of the state of the two-level system generally we 
have 

PnW = 1 -pfsW = l«f X (a| + (t))(+(t)|a) + |/3f x {a\ - (t))(-(t)|a) + 

x(a| + (t))(-(t)|a) X + /3V* x (a| - (t))(+(t)|a) x (75) 

and pf2(t) = |af x (a| + (t))(+(t)|6) + x (a| - (t))(-(t)|6) + 

x(a| + {t)){-{t)\b) X (<l>-(t)l<^>+W) + /3V* X (a| - (t))(+(t)|6) x ($+(t)|$_(t)). 

For the system initially prepared in one of the pointer states | ± (to)) (i.e. for a' = 
or (3' = 0) the above expressions can be simplified. For example, for the system initially 
prepared in the | + (to)) state generally we have 

PnW = M + and p^t) = {a\ + (t)).(6| + (t))* ; (76) 

while for the system initially prepared in the | — (to)) state we have 

pfiit) = \{a\ - m' and p^t) = {a\ - {t)).{b\ - (t))* . (77) 

Using equation (58) for the pointer states of our spin-boson model, the above 
equation reads 

PnW = cos2(^ + ^) and pf2(t) = -%in(<^ + ;^) for l^jj^^to)) = \ + (to)) 

pfiW = «m^(^-^) and pf,(t) = %in(^ - ;^) for |^^(to)) = | - (to))- (78) 
2 4\/n 2 2\/n 



The above expressions for pf2{t) basically mean that for the system initially prepared 
in one of its pointer states, the offdiagonal element of the reduced density matrix of the 
system should be a sinusoidal function with the frequency of Also, it must have 
successive zeros which are apart from each other by At = Airy/^/g. 

An examination of pf2(t) by plotting its more exact expression, given by equation 
(71), shows good agreement with the above result only as long as we have a very large 
average number of photons in the field. However, for a smaller average number of 
photons although we observe the oscillating behavior with the same frequency of 
we can clearly observe a decaying envelope which would destroy the offdiagonal element 
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of the reduced density matrix at large times (causing decoherence of the state of the 
system even when the system is initially prepared in one of its pointer states). Also, 
we observe that this decay is specifically more significant for a smaller average number 
of photons. Hence, we can guess that the difference between the prediction of equation 



(78) and what we expect from the more exact expression of equation (71), for the case 
that we have a smaller average number of photons, must be due to the fact that in 
calculating the pointer states of the system we assumed having a large average number 
of photons in the environment (so that we have a sharp distribution for the coherent 
state of the field and can assume y/n + 1 — y/n ^ ^^)- other words, we guess that 
the decoherence of the state of the central system when we start from one of the pointer 
states of the system must be due to having a limited number of photons in the field. 
In what follows our first goal is to make the appropriate corrections in equation 



(78) so that we can theoretically justify the decoherence of the state of the system when 



starting from one of the pointer states. Followed by that, we make corrections to the 



other elements of equation (75); and finally, we will use equation (75) together with the 
corrections which we make for having a limited average number of photons, in order to 
obtain a closed form for pf2(^)- ^^^^ s^^' after these corrections our closed form 

for the offdiagonal element of the reduced density matrix of the system will be in good 



agreement with the more exact but cumbersome expression of equation ( 71 ) which we 
obtained in this section for pf2(t). 

5.2. First order corrections due to having a finite average number of photons in the 
environment 

In section 3, while obtaining our pointer states of the system and the environment, 
we assumed having a large average number of photons in the environment; so that we 
could substitute \/n + 1 — -y/n by in our expressions. Now we consider the next 
order in the Taylor expansion of a/tT+T — ^/n about n, i.e. in 

t - VS) « E - + ...), (79) 

n=0 n=0 '^V'l' 

and make the appropriate corrections (due to having a finite average number of photons) 



in equation (78). In fact, by looking at equations (52) and (79) we notice that only at 
the limit of a large average number of photons, where Z^^o {y/n -\- 1 — y/n) ^ is 
a good approximation and there is no need to consider the next terms in our expansion 
for a/tT+T — i/n , the function G{t) will be independent of the states of the environment 
and pointer states can be realized for the system and the environment, which do not 
entangle with each other. 

We make corrections on pf2(t) of equation (78) by using the following substitution 



m our expressions 

oo 



n=0 
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Using equations (40) and (79) we have 



n=0 n=0 



e X 2^ ^ j — — = e "e x exp(?T, e^n^/a) (gi) 

n=0 ^• 



it' 3it' it' it' -it' 3it 

exp(r;, [e*"^''^ — 1]) e 4^/** = exp(n es"^/^ [esn^/^ — es^^]) e 



4V^ ■ 



which simphfies as 



u!_ . t' 



y |c„P e-**'(^^-^) ^ exp(2z n es^'^/^ sin(^— )) x e 4";%. (82) 
^0 8^3/2 

For an average number of photons large enough and times short enough for which 

t' 1 „,„ „„„ „„„„„„;™„ + „ „;„^ t' N U„ i' „„J U„ 1 i ii'^ 

n3/2 

can write 



< 1 we can approximate sin(^p72) by and em^^ by 1 + in the above 



equation. In other words, provided t goes to infinity slowly enough to have ^ 1 we 
Y: \cn\' e-^*'(v^-v^) ^ exp(2z n (1 + ^) x (^)) x e's^ or 



n=0 ^ 

00 , / 

+/2 /cni52 



Y |c„|2 e-**'{^"+i-^") ^ e-57ir e-*'^/32-^ (83) 

n=0 

So, to make the appropriate corrections in our expressions we should use the following 
substitution 



Making the above substitution in the expressions of equation (78) for pf2{t) we find 



pf2W = -^sin((^+^)^-^sin((^+^)e-*"/''^' for [^^(to)) = | + (t^)) and 



pUt) = sin(^ - ^) -> ^ - e-'^'^'^' for |^5(to)) = | - (to))- (85) 



One interesting aspect of the evolution of coherences given by the above equation 
is that for y9 = and = vr no matter whether the system initially is prepared in 
the I + (to)) state or the | — (to)) state, the evolution of pf2(t) is given by pf2(t) = 
=F| sin(^^) e~*'^/32n2 ^^^^j^ ^-^^ minus sign for = and the plus sign for = tt). Also, 
if V9 = I or y9 = ^, for both initial pointer states the evolution of |pf2(t)| is given by 
|Pi2(^)l = i'^°^(27H) e"*'^/^^"^ In general, for ip = mi/2 the evolution of |pf2(t)| will 
be the same for both initial pointer states. However, as we will see in the following 
paragraphs, this does not mean that for ip = mi/2 the evolution of |pf2(t)| becomes 
independent of the initial state of the system. 



In Figure 2 we have used equation (85) to plot the evolution of |pf2(t)| for the case 



that the system initially is prepared in the | + (to)) state. We also used the more exact 



expression, given by equation (71), to plot the same function. As we see, the correction 



(due to having a finite average number of photons in the environment) that we made on 
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Figure 2. Evolution of |/of2(^')l (wlierc t' = gt) for tlie case tliat the system initially 
is prepared in the | + (io)) state. Here we chose ip = tt/6 and n — 50. The curve 
represented by dashed lines is plotted by using the approximate expression which we 



obtained from our pointer states, given by equation (85 1. The other curve with solid 



lines is obtained from the more exact expression of equation ( 71 1 



our initial expression for pf2(t), nicely describes the decaying envelope in the evolution 
of coherences of the reduced system S, which is given by the factor e"*'^/^^"-^. 

pf2(t) has also been calculated for the simpler JC- model of quantum optics by Gea- 
Banacloche [S]. Gea-Banacloche showed that for the case that the two- level system is 
initially prepared in one of its pointer states the evolution of |pf2(^)| is simply given by 
|pf2(^)l ~ k e"*'^/^^"^ (which is valid for t' <^ n^/^); as opposed to what we have for our 
model, given by equation (85). 

As we observe, when the system is initially prepared in one of its pointer states the 
evolution of coherences of our model is characterized by an oscillating envelope, given 
by I I sin((y9+ ^^)|; while for the simpler JC-model of quantum optics we never observe 
such "decayo-sinusoidal" behavior in the evolution of coherences. 

Coming back to equation (75) for the general evolution of coherences of the 
reduced system, we also need to calculate the expressions for {a\ + (t)) {—(t)\b) and 
{a\ — {t)){+{t)\b), as well as the overlap between the pointer states of the environment 
($_(t)|$+(t)), for our generalized spin boson model. Using equation (58), we can 
evaluate (a| + (t)){—(t)\b) and (a| — (t)){-\-(t)\b) as follows 

t' . t' 



(a| + (t))(-(t)|6) 
{a\-{t)){+{t)\b) 



-zcos(- + 

2 4v^ 



cos — — 



^""^2 + vi^ 



sm 



4Vn 
t' . 



and 



(86) 



2 AVn 

However, here also we should make the appropriate correction due to having a finite 
average number of photons in the environment. Such correction can be made again by 
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using equation (84) to obtain 



(a| + {t)){-it)\b) = -.cos(^ + cos(^ - e-^'^"'^' and 

(a| - m^m) = ^sin(^ + sin(^ - e-*'^/«^^^ (87) 



Finally, we use the expressions for |$±(t)), given by equations (62 ) and ( [63j ), in order 
to calculate the overlap between the pointer states of the environment ($_(i(:)|$+(t)). 
We have 



|$±(t)) = e^¥(v/^+v^) |n) = 5: ^^^^ eT¥(v^+v^) |^). 

n=0 n=0 V^' 

As we discussed, for the coherent field with a large average number of photons we can 
use 

_ r- [n — n) (n — fiY , 



So, using the above relation and t' = gt, equation (88) becomes 

CO -n/2 -imp ^, (n + l-n) (n + l-n)^ r= {n-n) {n-nf . 

\^±{t)) ^ ' e^^^"^^^^ J^^sT^+V^+^TT-^;^) 1^) 



e-^/^eMA'V^+-^ - X E (90) 



2 ^4 ' 4v^ 8#/2J — vn! 

"^"P(4^Q"t2-4|]-4|7.» 
For the coherent field and within the approximation that we are using, as we 
discussed, J2n (n/n) ^ J2n |cnP = 1- Therefore, in the limit of n — )■ oo we can 
replace the expression [| — ^] of the above equation by |. So, using equation (40) we 
can simplify equation (90) to 

I*. W) « exp(T^|^ + « - 8^51"^ I " ^''Pfr^li - >• <»^' 

So now 

(<l>_(t)|$+(t)) ^ exp(-it'v/fi ["^ + — ^]) X ( z/ e2^'^"4kl| ^ g-syili-ik] (92) 

4 4n 8n"' 



Using the following formula from quantum optics [10] for the scalar product of the 
coherent states 

{u » = exp[-(|z/ f + \u\^)/2 + u'* u]. (93) 



equation (92) becomes 

(<l>_(t)|<l>+(t)) ^ exp(-it'v^ [7 + ^ - ^]) X exp( n {eW^t"?^) - 1}). (94) 

4 An on 

Therefore, 

|(<|._(t)|$+(t))P ^ exp(-4nsin2([^][^ - ^])). (95) 



Time- dependent pointer states of the generalized spin-boson model and etc. 



22 




Figure 3. Short time evolution of |pf2(i')l fo^' the case that the system initially is 
prepared in the lower state. Here we chose ip = n/Q and fi = 50. The curve represented 
by dashed lines is plotted by using the approximate expression which we obtained from 



our pointer states, given by equation (97). The other curve with solid lines is obtained 



from the more exact expression of equation ( 71 1 



For an average number of photons large enough and times short enough for which 
^ ^ 1 the above expression for the overlap between the pointer states of the 
environment reduces to 

.5 1 



\{^4m4tW^exp{-t" [ 



An' 



_ 25 ./2 
~ g 16 



(96) 



Now, using equations (75), (85), (87) and (94) we can obtain the following closed 
form for pf2(t), which is valid for t' <^ n^^'^ 
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In Figure 3 we used equation (97) to plot the short time evolution of |px2(^)l ^ 
case that the system initially is not prepared in one of its pointer states. We also used 
the more exact expression, given by equation (71), to plot the same function. As we see 
from this figure, equation (97) serves as a good approximation in closed form for the 
more exact relation, as long as we are not considering long times 

As we can see from equations (96) and (97), at sufficiently short times t' <^ ^Jft 



the decay of the first two terms is characterized by the decaying factor e"*'^/'^^""', while 
the decay of the other two terms is characterized by the much faster-decaying term 
due to the overlap between the pointer states of the environment ($_(t)|$_|_(t)) which 



^ For longer times, the approximate equation (971 shows internal oscillations in the evolution of |pf2(i)| 
which are misplaced compared to those of the plot which we obtain from the more exact expression of 
equation (71). However, the envelopes still do coincide with each other with great precession. 
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is proportional to the factor e~55*'^. This fact clearly shows why indeed we should 
generally expect a much slower decoherence of the state of the system when the system 
is initially prepared in one of its pointer states, compared to the case that the system 
initially is not in any of its pointer states. Also, equation (97) shows that for (p = mi/2 
and if the system initially is not prepared in one of its pointer states, the evolution of 
coherences of the central system would not be independent of the initial state of the 
system (just unlike the case that the system is initially prepared in one of its pointer 
states). However, since the last two terms of equation (97) vanish much faster than the 
first two terms, at larger times and for (p = mT/2 we expect the evolution of coherences 
of the central system to be independent of the initial state of the system. 

In this section we calculated the offdiagonal element of the reduced density matrix 
of the system in the basis of eigenstates of the operator. However, we could equally 
study the decoherence of the state of the central system in the basis of the | ± (to)) 
states. 

From equation (58) it is clear that for t' <^ the pointer states of the system 
almost are time- independent; and they can be approximated by |±(to))- So, in this limit 
the evolution of an arbitrary initial state of the central system [a' \ + (to)) + P' \ — (to))) 
in contact with an initial coherent field |z/) is approximately given by |'?/'*°*(t)) = 
a' I + (to))|$+(t)) + I - ito))\^-it)). Therefore, we would have 

ps{t) = |af I + (to))(+(to)| + 1/3? I - (to))(-(to)| 

I + (to))(-(to)| X (<l>_(t)|<l>+(t)) I - (to))(+(to)| X ($+(t)|$„(t)). (98) 

For this short range of times the evolution of the pointer states of the environment can 
be approximated by equation (91). So, in the | ± (to)) basis and for t' <^ we must 
have 

pf,(t) = a'/3-($_(t)|$4t)) 

3 3 1 , , =i£(5 

4 An 8n' 



a 13'* exp{-it'Vfi + jz: - ^]) x exp( n {e ^ ^* - 1}). (99) 



Finally, using equation (96) we find that for t' <^ A/n 

\pUt)\' ^ la'Pfexpi-t'^ - ^ la'Pf e'i^r (100) 

4 An 

Hence, in the basis of the | ± (to)) states the short-time decoherence of the state of 
the central system is characterized by the fast-decaying factor e~32*'^ when the system 
initially is not prepared in one of its pointer states; while in this basis the pointer states 
of the system almost do not decohere within short times. 



6. Summary and conclusions 

Considering a single-mode quantized field in exact resonance with the tunneling matrix 
element of the system, we obtained the time-evolution operator for our model. Using 
this time-evolution operator then we calculated the pointer states of the system and 
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the environment for the case that the environment is initially prepared in the coherent 
state with a large average number of photons. Most importantly, we observed that 
for our spin-boson model represented by the Hamiltonian of equation ([T]) the pointer 
states of the system turn out to become time-dependent; as opposed to the pointer 
states of a simplified spin-boson model (with Hg proportional to rather than a^) 
for which [Hs,Hmt] = 0. The simplified model has often been used in the context of 
quantum information and quantum computation to gain some insights regarding the 
decoherence of a single qubit flbl [T6l [T7] . However, in most of the practical situations 
different noncommutable perturbations may exist in the total Hamiltonian of a realistic 
system-environment model which would result in having time-dependent pointer states 
for the system Indeed, the authors believe that the fact that the pointer states of a 
system generally are time-dependent and may evolve with time has not been seriously 
acknowledged in the context of quantum computation and quantum information. In 
specific, in the context of quantum error correction [HI [19] it is often assumed that the 
premeasurement by the environment does not change the initial pointer states of the 
system. In other words, quantum "nondemolition" premeasurement by the environment 
is often assumed [ISl ES]; as is also assumed in Von Neumann's scheme of measurement 
[20] [1]. Also, in the context of Decoherence- Free-Subspaces (DFS) theory the models 
which often are studied either contain self-Hamiltonian for the system which commutes 
with the interaction between the system and the environment, or it is assumed that 
we are in the quantum measurement limit or in the quantum limit of decoherence 
flbl [2T1 122], [23]. However, all of these assumptions are in fact a big simplification of the 
problem; since, as we discussed in paper I, they completely exclude the possibility of 
having pointer states for the system which may depend on time jl]. 

Another interesting point in obtaining the pointer states of the system and the 
environment for our model was the realization of the fact that only in the limit of a 
large average number of photons can we have a set of (time-dependent) pointer states for 
the system. In other words, unless we have a sufficiently large average number of photons 
which can make a sharp distribution function for the state of the electromagnetic field, 
there is always some degree of entanglement between the states of the system and the 
environment (see equations (52) and (79)) and the pointer states of measurement cannot 
be realized at all. 



We also showed that at t = (2n+ l)7rv^/g (with n = 0, 1, 2, ...) the | ± (t)) pointer 
states of the system coincide with each other and hence, whatever is the initial state of 
the system, at these specific times the states of the system and the environment are not 
entangled with each other and the system can be represented by a well-defined state of 
its own. Using our pointer states, we also obtained a closed form for the offdiagonal 
element of the reduced density matrix of the system and studied the decoherence of 

* In the quantum measurement limit the interaction between the system and the environment is 
so strong as to dominate the evolution of the system H « i?int- Also in the quantum limit of 
decoherence the Hamiltonian for the system almost dominates the interaction between the system 
and the environment as well as the self-Hamiltonian of the environment H w H^. 
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the central system in our model. We showed that for the case that the system initially 
is prepared in one of its pointer states, the offdiagonal element of the reduced density 
matrix of the system will be a sinusoidal function with a slow decaying envelope which is 
characterized by a decay time proportional to n (through a decoherence factor calculated 
as e-g^*^/32n^~). ^iiile for the case that the system initially is not prepared in one of 
its initial pointer states, it will experience a fast decoherence within a time of order 
1/g. The "decayo-sinusoidal" evolution of coherences (figure 2) which we observe in 
our model and for the case that the system initially is prepared in one of its pointer 
states is a new form of decoherence which cannot be observed in the somewhat similar 
Jaynes-Cummings model of quantum optics [S]. 

It will be interesting to generalize this study to the case that the environment is 
not merely represented by a single-mode bosonic field; and consider some classes of 
spectral densities for the environment. Also, for the spin-boson model represented by 
the Hamiltonian of equation ([T]) at least in principle one should be able to obtain the 
pointer states of the system and the environment in some nonresonance regimes and for 
the single-mode quantized field. 

To further demonstrate the generality and usefulness of our method of obtaining 
pointer states, in another article [21] we will obtain the time-dependent pointer states 
of the system and the environment for the quantized atom-field model and in some 
nonresonance regimes. 



Appendix A 

In this section we introduce a measure for degree of entanglement between the states 
of the system and those of the environment, and by calculating this measure for the 



initial pointer states of our model (which are given by equation (55)) we show that our 
result for the pointer states of the system and the environment is valid over a length 
of time which is proportional to n, the average number of photons in the field; since as 
we will see, only within up to this range of times our pointer states of the system and 
the environment can stay separated and will not considerably entangle with the states 
of another subsystem. 

For the global state of the system and the environment, given by equation (??), we 

have 

= |A(t)) \a) + |B(t)) \b) 

= \A'{t)) [G{t)\a)] + \B{t)) lb); where \A'{t)) = (A-1) 

For our pointer states (given by | ± (t)) = Af± {G±{t)\a) + \b)}) we can define a degree 
of entanglement through the following relation 
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which also is equal to 



(A'(t)lB(t)) 



(A-3) 



|A'(t)P • 

The above function basically is the overlap between the vectors \A'{t)) and |B(t)) of 
equation (A-1); which is normalized to the unity, since for our pointer states of the 
system we have |A(t)) = G±{t)\'B{t)). From equation (A-1) it is clear that for perfect 
pointer states, where there is no entanglement between the states of the system and the 
environment, q{t) must always remain equal to the unity; all throughout the evolution 
of the system and the environment (i.e. \A'{t)) and |B(t)) must perfectly coincide with 
each other, in which case the states of the system and the environment in equation 
(A-1) will not entangle with each other; and we will have pointer states for the system 
which are given by equation (11)). Only in this case the states of the system and the 
environment in equation (A-1) will always stay separated and one can assign each of 
the two subsystems with well-defined states of their own. 

Our goal is to calculate our measure of entanglement q(t) for the pointer states of 
the system and the environment which we obtained for our model in this paper; and 
to study its evolution with time. For the global global state of the system and the 
environment, given by equation (A-1), one can calculate the reduced density matrix of 
the system as 

p^(t) = |a)(a|(A(t)|A(t)) + |6)(6|(B(t)|B(t)) + |a)(6|(B(t)|A(t)) + |6)(a|(A(t)|B(t)).(A-4) 
Therefore, 

pUt) = {B{t)\A{t)y, (A-5) 

where in the above equation pf2(^) the offdiagonal element of the reduced density 
matrix of the two- level system in the \a) and \b) basis. 

From equation (A-5) we can easily see that if the system initially is prepared in one 
of its initial pointer states; i.e. if lip^ito)) 



± (to)); then we have 
(A(t)|G±(t)B(t)) = Glit){B{t)\A{t)) = Gl{t)pUt) 
Therefore, 

(A(t)|G±(t)B(t)) Gl{t)pUt) 







\Mt)\ X 


B{t)\ 



\pf2{t)\ 



or 



^/i^X^l-pf,{t) 



(A-6) 



(A-7) 



(A-8) 



In calculating the above equations we must keep in mind that pfi(t) and pf2(t) must be 
calculated in the basis of the \a) and \b) basis states; and also they must be calculated 
for the case that the system initially is prepared in one of its initial pointer states, given 
by equation (55). 
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For our spin-boson model and for example for the case that the system initially is 



prepared in the | + (to)) state, from equations (78), (84) and (85) we had 



Pfi(t) = cos^(f + t4=) ^^"^^'^'^ \pUt)\ = hsini^ + ^)| e-*'V3-^ (A-9) 

2 4y/n 2 2\/n 

Therefore, 

|g(t)| = ^^J^^^^^^= ~ 1 if and only if t' < n (A-10) 
^pfi(t) X ^1-pUt) 

One can easily verify that for the case that the system initially is prepared in the | — (to)) 
state also, we would have the same result for the degree of entanglement between the 
states of the system and the environment. These results basically indicate that our 
result for the pointer states of the system and the environment is valid over a length 
of time which is proportional to n, the average number of photons in the field; since 
within times of the order of n/g the degree of entanglement, calculated for our pointer 
states, will stay close to the unity and our calculated pointer states will be immune to 
entanglement. 
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